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Abstract. Real-world networks, e. g. the social relations or world-wide-web graphs, 
exhibit both small-world and scale-free behaviour. We interpret lattice triangulations 
as planar graphs by identifying triangulation vertices with graph nodes and one- 
dimensional simplices with edges. Since these triangulations are ergodic with respect 
to a certain Pachner flip, applying different Monte-Carlo simulations enables us to 
calculate average properties of random triangulations, as well as canonical ensemble 
averages using an energy functional that is approximately the variance of the degree 
distribution. All considered triangulations have clustering coefficients comparable with 
real world graphs, for the canonical ensemble there are inverse temperatures with small 
shortest path length independent of system size. Tuning the inverse temperature to a 
quasi-critical value leads to an indication of scale-free behaviour for degrees k > 5. 
Using triangulations as a random graph model can improve the understanding of 
real-world networks, especially if the actual distance of the embedded nodes becomes 
important. 


Keywords: unimodular lattice triangulations, networks, maximal planar graphs 

PACS numbers: 64.60.aq, 02.10.Ox, 05.10.Ln 


Submitted to: New J. Phys. 



Unimodular lattice triangulations as small-world and scale-free random graphs 2 

1. Introduction 

Real world systems or networks that consist of many similar entities that are interacting 
can be described by graph theory (see DEI for reviews). Examples are the world-wide- 
web, where single websites are modelled as vertices and links between websites are 
modelled as edges; the network of social interactions, where vertices represent single 
humans and edges imply an existing friendship between two humans and the scientific 
co-author network, where authors of scientific papers are vertices and an edge between 
two authors is existent if both are co-authors of a common paper. 

Most of these real networks share three common properties 

(i) The clustering coefficient, which is basically the probability that two neighbours 
of a common node are connected by a graph edge, is high and independent of the 
size of the network. In real networks this means for example that there is a high 
probability that if two persons know a common person, they know each other, too; 
if two scientific authors are co-authors of a common third author, also these two 
are very likely co-authors. 

(ii) The length of the shortest path between two random vertices, measured by the 
number of edges the path consists of, is small even for huge graphs and scales 
like the logarithm of the system size for increasing system size. This behaviour is 
commonly known for the network of social interactions, where every two random 
persons of the world should know each other over less then 10 middlemen, or in 
the world wide web, where each website can be reached from every other website 
within a few clicks. 

(iii) The distribution of the number of edges incident with a vertex (its degree k) follows 
a power-law distribution ( P{k ) oc /c -7 ) with 2 < 7 < 3, which means that there are 
a lot of vertices with only a few neighbours, and only a minority highly connected 
vertices. For example in the co-author network there are few authors that worked 
together with a lot of different people, and there are a lot of authors that worked 
only with very few co-authors. 

The appearance of the first two properties is often denoted as small-world behaviour, 
networks with the third property are denoted as scale-free. 

To understand the structure and behaviour of real world networks different random 
graphs are used as model systems [3]. Three widely used models are the Erdos- 
Renyi random graph 0 El i, the Watts-Strogatz random graph [7] (and a slightly 
altered version denoted as Newman-Watts random graph 0 ) and the Barabasi-Albert 
random graph [9]. While all three models have a shortest path length scaling with 
the logarithm of the system size Emm], none exhibits both the other two properties of 
real-world networks: The Erdos-Renyi random graph has vanishing clustering coefficient 
and a binomial degree distribution; the Watts-Strogatz random graph exhibits a non¬ 
vanishing clustering coefficient for a special choice of parameters 0 , but also a degree 
distribution that agrees quantitatively with the Erdos-Renyi graph HH; the Barabasi- 
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Albert model creates a scale-free degree distribution with exponent 7 = 30 H21 US] , 
but leads also to a vanishing clustering coefficient for increasing system size [13] . Several 
more complicated models developed afterwards combine both small-world and scale-free 
behaviour mmm, there are also models that use graphs with vertices embedded in 
some geometric space that exhibit these two properties dHlEM. 

We present here a novel type of (embedded) random graphs by identifying 
triangulations of integer lattice point sets as graphs that shows a crossover from 
ordered, large-world to unordered, small-world and possible scale free behaviour. Such 
triangulations are tessellations of the convex hull of the point set into non-intersecting 
simplicial building blocks (triangles in two dimensions, tetrahedrons in three dimensions) 
EH- Triangulations are an important tool in physics for describing curved space (-times), 
in quantum geometry (e.g. in the framework of Causal Dynamical Triangulations 
m and in spin foams mv, they are also a major object of study in topology and 
geometry where one is for example interested in the number of distinct triangulations of 
a given topological manifold [23]. Triangulations are commonly used for describing foams 
[25], where often a special so-called Delaunay triangulation, which is the dual to the 
Voronoi tessellation, is used for foam construction. They can also be used for describing 
the topological properties of foams in terms of the neighbouring cells HH EH Eg. 
Additionally glass-like dynamics near regular configurations can be found in dual graphs 
of topological triangulations [29] . The dual graphs of disordered triangulations as 
described in this paper can be used for the construction of poly-disperse foams, which 
exhibit a broad range of cell sizes, and for considering e.g. their transport and diffusion 
problems. 

Triangulations have been used as random graph models before: Random Apollonian 
Networks |301 31, [32jJ (randomised dual graphs of Apollonian packed granular matter) 
are triangulations constructed using a procedure similar to preferential attachment and 
show both small-world and scale-free behaviour. Since each graph can be embedded into 
a closed surface with high enough genus, and triangulations are maximal planar graphs 
in the sense that each insertion of another edge will break planarity, triangulations of 
surfaces with arbitrary genus were studied in [33]. Canonical ensembles of triangulations 
of the sphere were used in [33j to consider a quench from the random triangulations to 
zero temperature. 

In contrast to the triangulations considered in [30] [31, [32] J33| [MJ before, which 
we denote as topological triangulations, we use embedded triangulations with vertices 
having fixed coordinates, precisely unimodular triangulations of two-dimensional 
integer point lattices. For topological triangulations only the topological degrees of 
freedoms (the way how vertices are connected) are important, contrary to embedded 
triangulations where additionally the actual coordinates of the vertices are fixed and 
specified. This causes complexities and numerical difficulties that do not have to be 
addressed in topological triangulations, e.g. to determine whether a Pachner flip leads 
to a valid triangulation. 

In this paper we measure the degree distribution, the clustering coefficient and the 
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Figure 1 . Examples of unimodular triangulations of a 8 x 8 lattice. The colour of 
the vertices corresponds to the number of incident vertices, the dashed border is fixed 
and will not be flipped, but is considered for the calculation of the energy. From left 
to right: Minimal energy triangulation (/3 —> +oo), random triangulation at infinite 
temperature (/3 = 0), example of a maximal energy configuration (/3 —► —oo). The 
colour code corresponds to the difference of the actual and the reference vertex degree 
as defined in ([!]). 


shortest path length of random lattice triangulation graphs using Metropolis Monte- 
Carlo simulations and find a high clustering coefficient and similarities to common 
network models. Introducing the notion of an energy of a triangulation that is well 
known in literature and corresponds to the variance of the degree distribution we apply 
the usual notion of statistical physics and examine the canonical ensemble averages of 
these observables for different values of the inverse temperature. For the numerical 
calculation of the expectation values we use the Wang-Landau algorithm for calculating 
the density of states. This makes it possible to calculate equilibrium properties for 
all temperatures, in contrast to Metropolis or Glauber dynamics used in |29, [33J 134] . 
where it is hard to access negative temperatures and low temperatures, but we loose 
the ability to consider dynamical behaviour like quenches or glass-transitions studied in 
the literature before. In the canonical ensemble we find in all considered observables 
a transition from an ordered large-world behaviour for positive temperatures to small- 
world and scale free behaviour for negative temperatures. 


2. Triangulations as random graphs 


A triangulation of a discrete point set A = {aq,... Xk}, Xi G l 2 is a tessellation of the 
convex hull of the point set into triangles <Ji E A 3 so that the interiors of two distinct 
triangles cq ^ <7j do not intersect, and so that every point A, is the corner of at least 
one triangle EB- Sometimes the latter property is not required, if all points are corners 
of a triangle the triangulation is then called full, if there are points that are not part of 
the triangulation it is called non-full. In this paper we consider only full triangulations 
of two-dimensional integer M x A-lattices 


A = 



m E Z M {0,1,..., M — 1}, n E Z 


N 
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which are unimodular, that is that all triangles have equal area 1/2. Examples of such 
triangulations can be found in £gure[lj 

There are analytical bounds [35] and numerical calculations [36] for the total 
number of triangulations on such integer lattices. Both the bounds and the numerical 
calculations show that number of triangulations scales exponentially with the system 
size of the underlying lattice, this extensivity makes it possible to apply tools from 
statistical physics to triangulations. A similar system was also considered in [37], 
where the convergence of a Metropolis-like Monte Carlo algorithm applied to lattice 
triangulations was analysed for different parameter choices. 

Taking the edges and the corners of triangulations in addition to triangles into 
account as simplices, a triangulation is also a simplicial complex, i.e. all faces of 
simplices are also contained in the simplicial complex. We interpret triangulations 
as graphs by identifying vertices as graph nodes and the triangulation edges as graph 
edges. Neglecting the boundary vertices triangulation graphs are maximal planar, i.e. 
no edge between an internal and another vertex can be inserted without violating the 
planarity of the graph. 

Since we are not interested in properties of certain triangulations but in averages 
over all triangulations of a given point set, we need a method that is able to construct 
all possible triangulations of a given point set. Therefore we use Pachner flips [38], 
elementary steps that map a triangulation into another one: Select an edge {xi,Xj} of 
the triangulation that connects two neighbouring triangles {x il Xj, Xk} and {xi,Xj,Xi}. 
If the quadrangle consisting of the two triangles is convex, replace the two triangles 
by {xi, Xk,xi} and {xj, Xk, x{\\ if the quadrangle has a concave or a flat angle, discard 
the flip as non-executable (see figure[2] for some examples). These flips are ergodic 
for triangulations of point sets in two dimensions [39] . that is that every triangulation 
of a given point set can be transformed into every other one by a finite number of 
Pachner flips, so we can use these Pachner moves for doing Markov chain Monte Carlo 
simulations. Notice that the introduced Pachner flip conserves the number of edges 
and the unimodularity of the triangulation. In principle there are more Pachner moves 
(inserting and removing vertices) beside flips, but these lead to non-full triangulations 
and thus are not considered here. In the context of graph theory on topological 
triangulations pn m the Pachner moved here is called Tj, and the two additional 
moves not used here are called T 2 . Note that there is a fundamental difference between 
the diagonal edge flips in topological and in lattice triangulations: Each former flip can 
always been executed, whereas latter ones can only be executed if the two adjacent 
triangles form a convex quadrangle. This is due to the fixation of the vertex coordinates 
in lattice triangulations, while for topological triangulations no coordinates of the 
vertices are specified and the notions of convex or concave do not exist. 

To categorise triangulations of a point set according to their order and disorder, we 
define for a triangulation T the energy function 

E(T) := £ W T| - y™) 2 

dGA 


( 1 ) 
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Figure 2. The two different Pachner moves in two-dimensional lattice triangulations. 
(Left) the diagonal edge flip that replaces one diagonal of a quadrilateral with the other 
diagonal. On the lattice this steps leads to a valid triangulation if the two triangles 
incident with the edge to be flipped constitute a parallelogram. (Right) the removal 
(up to down) and the insertion step, which are not used in our system, since we consider 
only full triangulations that are ergodic with respect to the diagonal edge flip. 


taking low values for ordered triangulations and high values for unordered triangulations, 
ffere k\, ' is the number of edges that are incident with the vertex v in triangulation T 
and 7o is a reference triangulation that will be the ground-state of the energy function 
(in this paper we choose the maximal ordered triangulation displayed in figure!!] as the 
reference triangulation). The additional term kf 0 with To being the maximal ordered 
triangulation can be seen as an implementation of fixed boundary conditions as displayed 
in figure]!] Periodic boundary conditions, which yield toroidal topology, cannot be used 
for lattice triangulations, because on closed surfaces without boundaries there are several 
possibilities to connect two vertices with an edge (e.g. with different winding number) 
and therewith it depends on the choice of the connection whether the considered object 
is a valid triangulation. In contrast, for the considered topology, the edge between two 
vertices is always the line segment. 

Since the number of edges of full lattice triangulations is constant, sums over 
linear vertex degrees vanish and the chosen energy function (]!]) is the simplest possible 
polynomial energy function in the number of vertices, edges and triangles. It can be 
related to the square of the local curvature usually used in dynamical triangulations 
[22] (which is basically the deficit angle given by the number of triangles around a 
vertex minus six). Similar energy functions were already applied to graphs in na ; in 
[37] the total length of all edges was used as energy function in lattice triangulations, 
which qualitatively agrees with our choice since high energy leads to long edges in 
the triangulations; previous works [29], [33], 33] considering topological triangulations as 
graphs also use the energy function ([lj with the mean vertex degree ( k) instead of kf£°\ 

Interpreting the triangulation vertices as graph nodes and the 1-simplices as graph 
edges, we compare the graph properties of the triangulations with Erdos-Renyi jHEllS], 
Newman-Watts [8] and Barabasi-Albert [9j random graphs. In all three cases we choose 
the model parameters so that the number n of the random graph nodes matches the 
number MN of the lattice triangulation vertices and the average number of random 
graph edges matches the number e = 3 MN — 2(M + N) + 1 of the lattice triangulation 
edges. 

In the Erdos-Renyi random graph there are n vertices, each pair of vertices is 
connected by an edge with probability p (U El E]- They have a small world behaviour 
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for the shortest path length, but vanishing clustering coefficient and a binomial degree 
distribution. To compare the Erdos-Renyi random graph with triangulations we choose 
its parameters to be n = MN and p = 2 e/n{n — 1) —> 6/n. 

The Newman-Watts random graph [8] is a modification of the Watts-Strogatz graph 
[7]. Starting with a regular graph of n vertices and connections to the next L neighbours, 
for each present edge an additional edge is inserted between two random vertices with 
probability q. This model can be seen as the periodic regular graph superimposed with a 
random graph. The random rewiring leads to shortcuts and a short average path length, 
the basic regular graph leads to a high clustering coefficient. To obtain a comparable 
random graph we use the parameters n = MN, L = 4 and q = — 1 + e/2 MN 

The Barabasi-Albert random graph [9] is constructed as following: Start with m 
isolated vertices and iteratively insert t vertices, each with m edges to already present 
vertices, such that the probability for connecting to a present vertex is proportional 
to its degree (preferential attachment). This random graph has vanishing clustering 
coefficient, but a power-law degree distribution and a small-world shortest path length. 
We choose the model parameters m — 3 and t = MN — 3 to compare with M x N 
lattice triangulations. 

3. Random triangulations 

In this section we calculate the scaling behaviour of the vertex degree distribution, 
the mean energy, the average clustering coefficient and the average shortest path 
length in terms of the system size of random triangulations and compare it with 
random graphs. With random triangulations we denote the ensemble with each possible 
triangulation having the same constant weight and contributes equally in the calculation 
of expectation values. The results are especially interesting, because due to the constant 
ensemble weights (that are in particular independent of the triangulation energies) they 
do not depend on the actual choice of the energy function. 

We calculate these observables in the random ensemble using the Metropolis 
algorithm [41] taking into account at least 1000 different random triangulations per 
considered system size. The autocorrelation time of all considered observables is smaller 
than 10MN Pachner moves, to avoid errors due to the autocorrelation, we execute 
1000M7V flips before the Erst and between two successive measurements to calculate 
the observable averages. 

Random triangulations were considered up to a system size of 64 x 64, which 
seems to be small if comparing with topological triangulations [34, [33]. For lattice 
triangulations one has to check the convexity of quadrangles to decide whether a step is 
permissible, which increases the computation time needed for one Metropolis step and 
which decreases the step acceptance ratio, because some steps have to be rejected. 

In this paper we restrict to quadratic integer lattices for clarity reasons, actually also 
simulations for non-quadratic lattices were performed. Except of the average shortest 
path length (which grows with |M — 7V|) all observables only depend on the actual 


JJnimodular lattice triangulations as small-world and scale-free random graphs 



vertex degree k system size MN 



Figure 3. Scaling behaviour of random lattice triangulations and comparable random 
graphs, a) Degree distribution for random triangulations of 8 x 8 (■) and 32 x 32 (♦) 
lattices compared to the average degree distribution of a Newman-Watts (□) and a 
Barabasi-Albert («) random graph with n = 32 2 vertices, b) Specific mean energy, 
variance of the degree distribution, c) clustering coefficient and d) shortest path 
length for maximal ordered (•*■) and random (•) triangulations, the Erdos-Renyi (o) 
and the other random graphs in terms of the number of vertices MN. For the random 
triangulations and the Newman-Watts graph the lines were determined by a power-law 
fit, all other lines are analytical results. The asymptotic behaviour is listed in Tab. [I] 


system size M ■ N for M ■ N <1 10 2 and M, N 4 and not on the two linear sizes. 

For comparing the random triangulation averages with the results in the Erdos- 
Renyi, the Newman-Watts and the Barabasi-Albert random graphs of the same size 
that are not known analytically, we use averages over 500 randomly generated instances 
of the respective random graphs. The generation of these random graphs was done using 
the NetworkX framework m- 

3.1. Vertex degree distribution 

The vertex degree distribution of the considered random graphs is known analytically: 
for the Erdos-Renyi and Newman-Watts it follows a binomial distribution, for the 
Barabasi-Albert random graph it is a power-law distribution [121113]. 

In figure^) the degree distribution of random triangulations is displayed for 
different lattice sizes and compared with the different random graphs. In contrast to 
the random graphs the values k — 0 and k — 1 are not encountered for triangulations 
since these vertex degrees are forbidden for the graph to form a valid triangulation. The 
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degree k — 2 can only be realised on the boundary of the triangulation, so the probability 
decreases for increasing lattice size and decreasing importance of the boundary. 

The calculations for our random triangulations agree for k > 4 with the results for 
the degree distribution of a topological triangulation of a torus in [SI]- In contrast to 
our results, where P(k) is peaked at k — 4, in [33] a monotonically decreasing P(k ) with 
maximal value at k — 3 was found. This is because the vertex degree k — 3 is difficult to 
realise due to the non-general coordinates of the vertices in lattice triangulations (there 
are many collinear points), so even by constructing one cannot realise more than one 
out of four vertices having degree k = 3. In topological triangulations one can realise 
Apollonian-like networks [30 , . 311132] , where every second vertex can have degree k — 3. 

To compare a Newman-Watts random graph we used a ring with each vertex 
connected to its L = 4 nearest neighbours and added edges, so for this graph the degrees 
k < 4 are not present. As a result one can see that the degree distribution of random 
triangulations is comparable with the Erdos-Renyi and Newman-Watts random graph 
for k > 5, whereas the Barabasi-Albert model shows a scale free power law behaviour 
oc k~ 3 . 


3.2. Mean energy 


To quantify the disorder of a triangulation we use an energy function which is, up 
to boundary terms, the variance of the degree distribution multiplied with the lattice 
size. In figure]^)) we compare the specific average energy of random triangulations with 
the degree variance of the random graphs, which can be calculated from the vertex 
degree distributions. For the random triangulations a fit to the numerical calculated 
data yields a convergence of the specific average energy to 4.62 ± 0.05, the Erdos- 
Renyi and Newman-Watts random graphs show a similar behaviour (see Tab. [I] for the 
detailed values). For the Barabasi-Albert random graph the second moment of the 
degree distribution diverges. 

For random topological triangulations the energy per vertex can be calculated 
analytically to be E/n = ({k) — 3) 2 = 9 for (k) = 6 in topological triangulations [33], 
which agrees with our finding that E/MN converges for MN —» oo. The absolute value 
of the specific energy in lattice triangulations is smaller, because the degree distribution 


has a peak located nearer at the mean value as discussed in Sec. 3.1 


3.3. Clustering coefficient 

The clustering coefficient Ci = 2K/ki(ki~ 1) of a vertex i with degree ki is the ratio of the 
number of connections K between the ki neighbours and number of possible connections. 
For a non-boundary vertex of a triangulation holds K = ki and Ci = 2/{ki — 1), for a 
boundary vertex K — k t — 1 and Ci = 2/fcj. The clustering coefficient C of the whole 
graph is the average of the vertex clustering coefficients. 

The clustering coefficient for the Erdos-Renyi random graph equals the edge 
connection probability p (~ 6/ MN to ensure equal edge number), for the Barabasi- 
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Albert model it can be calculated using a mean-field approach ra- In both cases the 
clustering coefficient vanishes for increasing lattice sizes for our choice of parameters. 

The numerical results for the random triangulations and the Newman-Watts graph 
can be seen in figure^), both are converging to a constant value for increasing lattice 
sizes. For random triangulations the limit is 0.4859 ± 0.0005, which is higher than 
for all considered random graphs (see Tab. [I] for the detailed values). Intuitively this 
means that almost every second possible edge between neighbours of a common vertex 
is present in the random triangulation. 

For topological triangulations the clustering coefficient is C ~ 0.6 [33]. This is 
higher than for our lattice triangulations, mainly due to the fact that the degree k = 3 
is much more probable for topological than for random triangulations. 

3.4. Shortest path length 

As for the clustering coefficient the shortest path length is analytically known for both 
the Erdos-Renyi and the Barabasi-Albert random graphs [TD] and has to be calculated 
numerically for the Newman-Watts graph. For all random graphs the shortest path 
length shows small-world behaviour and increases approximately with the logarithm of 
the vertex number. 

These values are compared with the shortest path length of random triangulations 
in figure[3ji). I 11 contrast to the random graphs, for random triangulations there is a 
power law behaviour of the shortest path length oc (Mjv) 0 - 400 ^ 0 - 001 . Also for random 
topological triangulations a power law scaling of the average shortest path length can 
be found [53]. 

4. Canonical triangulations 

In this section we consider canonical averages of the degree distribution, the mean 
energy, the clustering coefficient and the shortest path length of lattice triangulations 
in terms of the inverse temperature with respect to the energy defined in 0- These 
calculations extend the results obtained in the previous section which correspond to 
the special case of vanishing inverse temperature j3 = 0. The numerical simulations 
are not restricted to the usually considered case /3 > 0, but can be extended to the 
case of negative inverse temperatures (3 < 0. These negative temperatures can also be 
interpreted as positive inverse temperatures with a negative coupling constant included 
in the definition ([!]) of the energy, which makes the triangulations with the most disorder 
the ground state of the energy function. 

It is difficult to use Metropolis Monte-Carlo simulations BU to perform canonical 
averages for triangulations, especially if one wants to explore the regime of negative 
temperatures/coupling. There exist triangulations that are local minima in the energy 
landscape [36], and due to the small Metropolis acceptance probabilities (oc exp(— f3AE), 
where A E is the energy difference induced by the flip) the algorithm gets stuck in one 
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Table 1. Scaling and functional dependence obtained by a least-square fit of the 
mean energy (or the degree distribution variance), the clustering coefficient and the 
shortest path length with the system size for ground state (maximal ordered) and 
random triangulations as well as for comparable size Erdos-Renyi, Newman-Watts 
and Barabasi-Albert graphs. If there are no error bars, the results are taken from 
analytical calculations. 

quantity graph scaling behaviour 



ground state triang. 

0 


random triangs. 

(4.62 ± 0.02) - (10.80 ± 0.043) • (MN)-° A7±0 '° 2 

Var (k) 

Erdos-Renyi 

6 - 4 • ( MN ) 0 ' 5 


Newman-Watts 

2 - 4 ■ (MN) 0 ' 5 


Barabasi- Albert 

—> oo 


ground state triang. 

0.4 


random triangs. 

(0.4849 ± 0.0003) - (0.78 ± 0.02) • (M N)-°' 73±0 ' m 

(C) 

Erdos-Renyi 

6 • (MN)- 1 


Newman-Watts 

(0.241 ± 0.001) - (1.18 ± 0.02) • (mN)~°' 557±om7 


Barabasi- Albert 

0.678 • (MN)- 1 log(MN) 2 


ground state triang. 

(0.5677 ± 0.0002) ■ (. MN ) 0A9978±0 '° 0005 


random triangs. 

(0.673 ± 0.006) • (MN) 0A09±0 ' 0m 

<*> 

Erdos-Renyi 

\og(MN)/6 


Newman-Watts 

(0.063 ± 0.001) • log [(0.595 ± 0.002)MTV] 


Barabasi- Albert 

log(MN)/log(log(MN)) 


of these local minima. Other states outside of the local minimum that contribute to the 
ensemble average are reached only after many steps or not at all, so the autocorrelation 
times and therewith the simulation errors become large, or the possibility arises that 
the system is not even computationally ergodic anymore. This problem was also treated 
analytically in [37], where the mixing time (which is related to the autocorrelation time) 
of Glauber dynamics on lattice triangulations is shown to scale exponentially with the 
system size for a small enough f3 < 0. 

Using a parallel tempering approach [03], which is basically the parallel calculation 
of multiple Metropolis simulations at different inverse temperatures with the possibility 
of interchanging the inverse temperatures, can help to overcome the problem of these 
local minima in the energy landscape. But parallel tempering is known to fail in 
situations with a large free energy barrier, e.g. in first order phase transitions. 

Flat histogram algorithms as the Wang-Landau algorithm (H 05] are used to 
calculate the density of states (DOS) g(E) (which is the normalised number of states 
with energy E ) of systems and can help to overcome both the local minima and the 
free energy barrier problem. This algorithm samples the single states according to their 
inverse DOS g(E)- 1 based on an initial estimation of g(E) and simultaneously improves 
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Figure 4. Temperature dependence of graph properties of lattice triangulations, 
a) Canonical averaged degree distribution of 8 x 8 triangulations for different inverse 
temperatures /3. The dashed line is a power law oc k~ 3 . b) Mean energy in units of the 
energy of a random triangulation (with inset for positive temperatures), c) clustering 
coefficient and d) shortest path length in terms of the inverse temperature /3 for 
triangulations of different lattice sizes. The degree distribution and the lines in b)- 
d) were calculated using a Wang-Landau simulation, the data points in c)-d) are taken 
from a parallel tempering simulation. 


the estimation until it converges to the actual DOS of the system. The two main 
advantages of this algorithm is that with the knowledge of the DOS the observables can 
be calculated for all inverse temperatures f3 based on only one simulation, whereas for 
the Metropolis algorithm a new simulation for each f3 has to be used. The second one 
is that the algorithm does not get stuck anymore in local minima and the problems of 
large autocorrelation times described before and in [37] do not occur. Flat histogram 
sampling techniques were already applied to two-dimensional lattice triangulations [3 6] , 
but are limited to relatively small system sizes. 

In figure [5] one can see a comparison of Metropolis sampling and a sampling based 
on the DOS calculated with Wang-Landau in terms of the acceptance ratios and the 
autocorrelation time of the energy observable. For negative temperatures the Metropolis 
algorithm basically gets stuck in high-energy states, which leads to low acceptance ratios 
and long autocorrelation times. Using the Wang-Landau algorithm this problem does 
not occur, because steps are weighted according to their entropy difference, and not 
their energy difference. 

The DOS of two-dimensional lattice triangulations can only be calculated up to 
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Figure 5. Acceptance ratio of a Monte Carlo step (a) and autocorrelation time of the 
energy observable (b) in terms of the inverse temperature for the Metropolis algorithm 
for different lattice sizes (points). The solid line is the temperature independent value 
obtained by a multicanonical simulation. Note that for the simulations based on the 
density of states the acceptance ratio is lower and the autocorrelation time is higher 
than in the Metropolis case for certain temperatures ranges, but the Wang-Landau 
algorithm has on the one hand the advantage that only one simulation has to be 
performed for the whole temperature range, on the other hand it makes negative 
temperatures accessible at all. 


11 x 11 triangulations for all energies, introducing a energy cutoff makes it possible to 
calculate the density of states up to 25 x 25 triangulations that can be used only for 
positive inverse temperatures (d > 0. For bigger lattice sizes the Wang-Landau algorithm 
does not converge anymore in a reasonable amount of time, since on the one hand the 
entropy between neighbouring states differs by orders of magnitude, on the other hand 
the high energy states are only connected among themselves by pathes involving low 
energy states, which makes flat histogram sampling difficult (see [SB] for a detailed 
discussion). 

4-1. Degree distribution 

The canonical averaged distribution P{k ) of the vertex degrees of 8 x 8 lattice 
triangulations is plotted for different temperatures in figure[4|r). For the positive inverse 
temperature fd = 0.2 one finds a similar behaviour as for the random triangulations 
fd = 0 with the maximum of the degree distribution shifted towards degrees k = 5 
and k = 6, because this degrees are preferred at low temperatures due to the energy 
function ([Tj). For negative inverse temperatures fd = —0.1 there is an additional peak 
at vertex degree k & 40, because for this inverse temperature the triangulation with 
maximal energy contributes most to the ensemble average. An interesting behaviour 
can be found for the inverse temperature fd = —0.04 (here the probability distribution 
of the triangulation energies has two peaks, which is a hint for a phase transition). For 
this inverse temperature the degree distribution behaves similar to a power law k~ 3 
for vertex degrees in the range 5 < k < 25, which matches the degree distribution 
of the Barabasi-Albert model qualitatively. One has to admit that due to the small 
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system size the range of vertex degrees found spans less than one magnitude and that 
there are deviations from a pure power law behaviour even in this range. Nevertheless 
the qualitative change of the degree distribution compared to random triangulations is 
remarkable. Additionally we suspect that for bigger system sizes there is a clearer power 
law behaviour for larger ranges of degrees. 

Similar results for the canonical averaged triangulation degree distribution can be 
found in [34] for positive temperatures. 

4-2. Mean energy, clustering coefficient and shortest path length 

To compare the temperature dependence of the mean energy for different lattice sizes, 
we consider the expectation value of the mean energy e = (E) c (f3) /(E) rnd in terms of 
the energy of a random triangulation. The results of the numerical calculations for 
the relative mean energy, the clustering coefficient and the shortest path length can be 
found in figure]^)) - d). 

The temperature dependence of all considered observables shows the characteristic 
behaviour of a first-order phase transition, with quasi-critical inverse temperature 
/3 C —> 0 for increasing system size, but this is not a real phase transition. For the 
maximal energy of a quadratic N x N lattice triangulation there is a lower bound that 
scales with N 4 [36], so that the maximal specific energy scales at least with N 2 . For 
(3 < 0 (which is equivalent to using negative coupling) these maximal energy states 
are the ground states, which are then not bounded by below for the limit of infinite 
system size. For negative temperatures only finite systems can be considered and and 
the thermodynamic limit of infinite system size cannot be obtained, so the behaviour 
found at fd —>• 0 is then no actual phase transition in sense of statistical physics. 

For negative inverse temperatures one finds a high clustering coefficient 
(approximately 0.7 independent of the system size), and small average path length 
between 2 and 3 independent of the system size. The latter can be understood in terms 
of the graph theoretical diameter of the highest energy triangulations, which can be 
shown to be 4 independent of the lattice size. So for a small enough negative inverse 
temperature one can find a small-world behaviour for the triangulation graphs. 

For positive temperatures the average energy per vertex and the clustering 
coefficient can be approximated analytically. The first possibility is to use the mean- 
field approach suggested in [33], where the topological and entropic properties of the 
triangulations are neglected. This results in 

, cw 5e^/6 + 2/6 
MN ~ 2 + eH ' ' ~ 2e-P + 1 

The second possibility is to use the known degeneracy Q(E = 4) and Q(E = 0) of the 
ground state and the first excited state [3H] to calculate the expectation values neglecting 
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inverse temperature /3 inverse temperature /3 


Figure 6. Analytical approximations for the specific energy (a) and the clustering 
coefficient (b) in two different regimes compared with the actual data obtained 
numerically. The mean-field approach (red, dashed line), which is valid for an 
intermediate regime of positive inverse temperatures, neglects the entropic properties 
of the triangulations, the two-state approximation (blue, dotted line), which is valid 
for /3 —»• oo, just considers the ground state and the first excited state of the lattice 
triangulations. 


all other states, resulting in 

(E) _ n(E = 4) 4 _ 4 MN 

~MN ~ MN e 4/3 + Q(E = 4) ~ 2e ^ + M 2 N 2 

1 + (15MJV)->n(B = 4) e -« 

' ' » l + Sl(E = i)e-V 

Both approximations are compared with the numerical results in hgure[6j For both the 
specific energy and the clustering coefficient the mean held approximation is correct for 
0 < (3 < 1.5, the two-niveau low-temperature approximations is correct for fd > 3. 

5. Conclusion 

In this paper we proposed a new model for real-world graphs or networks by interpreting 
unimodular triangulations of two-dimensional integer lattice as graphs. Considering 
averages of random triangulations we found that they show a higher clustering coefficient 
than common random graph models, but also a power law growing shortest path length 
in terms of the system size, where one would expect a logarithmic scaling for small 
world behaviour. The degree distribution behaves similar as in the Erdos-Renyi and 
Newman-Watts random graph. 

Introducing an energy function that measures the order and disorder of a 
triangulation, canonical averages of graph observables in triangulations can be calculated 
using the Wang-Landau algorithm. This Monte Carlo method made it possible for 
the first time to calculate expectation values for the whole temperature range, albeit 
only for small system sizes. Using the Wang-Landau algorithm it is possible to 
access equilibrium properties of both high energy triangulations, which are poorly 
connected with other similar triangulations, and low energy triangulations, where 
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dynamically a glass-like phase was found in the literature before, both phenomena 
making it difficult to use the Metropolis algorithm. Considering the temperature 
dependence of the clustering coefficient, the average shortest path length and the degree 
distribution we found small-world behaviour for negative temperatures (3 < (3 C below a 
negative, quasi-critical temperature and hints for scale-free behaviour for temperatures 
near the quasi-critical temperature. All considered observables show a cross-over 
behaviour going from negative temperatures, where the partition function dominated 
by disordered triangulations, to positive temperatures, where it is dominated by ordered 
triangulations. 

In contrast to the topological the lattice or in general embedded triangulations 
are useful if one deals with networks and graphs where the actual coordinates of 
the vertices become important. Despite the introduction of vertex coordinates which 
induce non-executable flips, our results show qualitatively agreement with the graph 
properties of topological triangulations of surfaces with suitable genus that were 
estimated analytically and calculated numerically in the literature. 

An extension to the grandcanonical ensemble, where the number of vertices is not 
fixed any more, is possible by considering also insertion and removal Pachner moves. 
Additionally one can apply the methods described in this paper also to triangulations 
of arbitrary point sets in two and more dimensions. 

Acknowledgements 

The authors thank J. F. Knauf, R. Jonsson and A. Kraynik for fruitful discussions. This 
work is supported by EFI Quantum Geometry and the Elite Network of Bavaria. 

[1] Newman M 2010 Networks. An introduction (Oxford University Press) 

[2] Costa L d F, Oliveira O N, Travieso G, Rodrigues F A, Villas Boas P R, Antiqueira L, Viana M P 

and Correa Rocha L E 2011 Adv. Phys. 60 329-412 

[3] Albert R and Barabasi A L 2002 Rev. Mod. Phys. 74 47-97 

[4] Erdos P and Renyi A 1959 Publ. Math. Debrecen 6 290-297 

[5] Erdos P and Renyi A 1960 On the evolution of random graphs Publication of the Mathematical 

Institute of the Hungarian Academy of Sciences pp 17-61 

[6] Erdos P and Renyi A 1964 Acta Math. Acad. Sci. H. 12 261-267 

[7] Watts D J and Strogatz S H 1998 Nature 393 440-442 

[8] Newman M and Watts D 1999 Phys. Lett. A 263 341 - 346 

[9] Barabsi A L and Albert R 1999 Science 286 509-512 

[10] Fronczak A, Fronczak P and Holyst J A 2004 Phys. Rev. E 70 056110 

[11] Barrat A and Weigt M 2000 Eur. Phys. J. B 13 547-560 

[12] Dorogovtsev S N, Mendes J F F and Samukhin A N 2000 Phys. Rev. Lett. 85 4633-4636 

[13] Krapivsky P L, Redner S and Leyvraz F 2000 Phys. Rev. Lett. 85 4629-4632 

[14] Fronczak A, Fronczak P and Holyst J A 2003 Phys. Rev. E 68 046126 

[15] Holme P and Kim B J 2002 Phys. Rev. E 65 026107 

[16] Klemm K and Eguiluz V M 2002 Phys. Rev. E 65 036123 

[17] Szabo G, Alava M and Kertesz J 2003 Phys. Rev. E 67 056102 

[18] Rozenfeld A F, Cohen R, ben Avraham D and Havlin S 2002 Phys. Rev. Lett. 89 218701 

[19] ben Avraham D, Rozenfeld A F, Cohen R and Havlin S 2003 Phys. A 330 107 - 116 

[20] Xu X J, Wu Z X and Chen G 2007 Phys. A 377 125 - 130 



Unimodular lattice triangulations as small-world and scale-free random graphs 


17 


[21] De Loera J A, Rambau J and Santos F 2010 Triangulations. Structures for Algorithms and 

Applications (Springer) 

[22] Arnbjorn J, Jurkiewicz J and Loll R 2005 Phys. Rev. D 72 064014 

[23] Rovelli C 2007 Quantum Gravity Cambridge Monographs on Mathematical Physics (Cambridge 

University Press) 

[24] Sulanke T and Lutz F H 2009 Eur. J. Comb. 30 1965-1979 

[25] Sullivan J M 1999 The geometry of bubbles and foams Foams and Emulsions (NATO Science 

Series E vol 354) ed Sadoc J and Rivier N (Springer) pp 379-402 

[26] Oguey C, Rivier N and Aste T 2003 Eur. Phys. J. B 33 447-455 

[27] Aste T and Rivier N 1995 J. Phys. A 28 1381 

[28] Dubertret B, Rivier N and Peshkin M A 1998 J. Phys. A. 31 879 

[29] Aste T and Sherrington D 1999 J. Phys. A 32 70497056 

[30] Andrade J S, Herrmann H J, Andrade R F S and da Silva L R 2005 Phys. Rev. Lett. 94 018702 

[31] Zhou T, Yan G and Wang B H 2005 Phys. Rev. E 71 046141 

[32] Song W M, Di Matteo T and Aste T 2012 Phys. Rev. E 85 046115 

[33] Aste T, Gramatica R and Di Matteo T 2012 Phys. Rev. E 86 036109 

[34] Kownacki J P 2004 Eur. Phys. J. B 38 485-494 

[35] Kaibel V and Ziegler G M 2003 Lond. Math. Soc. Led. Note Ser. 307 277-308 

[36] Knauf J F, Kruger B and Mecke K 2014 Entropy of unimodular lattice triangulations 

arxiv: 1412.2541 

[37] Caputo P, Martinelli F, Sinclair A and Stauffer A 2013 Random lattice triangulations: Structure 

and algorithms Proceedings of the Forty-fifth Annual A CM Symposium on Theory of Computing 
pp 615-624 

[38] Pachner U 1986 Abh. Math. Sem. Univ. Hamburg 57 69-85 

[39] Lawson C L 1972 Discrete Math. 3 365 - 372 

[40] Farkas I, Derenyi I, Palla G and Vicsek T 2004 Equilibrium statistical mechanics of network 

structures Complex Networks ( Lecture Notes in Physics vol 650) ed Ben-Naim E, Frauenfelder 
H and Toroczkai Z (Springer) pp 163-187 

[41] Metropolis N, Rosenbluth A W, Rosenbluth M N, Teller A H and Teller E 1953 J. Chem. Phys. 

21 1087 

[42] Hagberg A A, Schult D A and Swart P J 2008 Exploring network structure, dynamics, and function 

using networkx Proceedings of the 7th Python in Science Conference ed Varoquaux G, Vaught 
T and Millman J (Pasadena, CA USA) pp 11 15 

[43] Earl D J and Deem M W 2005 Phys. Chem. Chem. Phys. 7 3910-3916 

[44] Wang F and Landau D P 2001 Phys. Rev. Lett. 86(10) 2050-2053 

[45] Wang F and Landau D P 2001 Phys. Rev. E 64(5) 056101 



